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ABSTRACT 

The magnetic field topology in the surrounding of neutron stars is one of the key 
questions in pulsar magnetospheric physics. A very extensive literature exists about 
the assumption of a dipolar magnetic field but very little progress has been made 
in attempts to include multipolar components in a self-consistent way. In this paper, 
we study the effect of multipolar electromagnetic fields anchored in the star. We give 
exact analytical solutions in closed form for any order l and apply them to the retarded 
point quadrupole (l = 2), hexapole (Z = 3) and octopole (Z = 4), a generalization of 
the retarded point dipole (Z = 1). We also compare the Poynting flux from each 
multipole and show that the spin down luminosity depends on the ratio R/r l, R 
being the neutron star radius and tl the light-cylinder radius. Therefore the braking 
index also depends on R/r l. As such multipole fields possess very different topology, 
most importantly smaller length scales compared to the dipolar field, especially close 
to the neutron star, we investigate the deformation of the polar caps induced by 
these multipolar fields. Such fields could have a strong impact on the interpretation 
of the pulsed radio emission suspected to emanate from these polar caps as well as 
on the inferred geometry deduced from the high-energy light-curve fitting and on the 
magnetic field strength. Discrepancies between the two-pole caustic model and our new 
multipole-caustic model are emphasized with the quadrupole field. To this respect, we 
demonstrate that working with only a dipole field can be very misleading. 

Key words: magnetic fields - methods: analytical - stars: neutron - stars: rotation - 
pulsars: general 


1 INTRODUCTION 


Rotating magnetized neutron stars are at the base of several compact object classes such as pulsars, magnetars, isolated 
neutron stars and X-ray binaries among others. Whereas the rotation period P and its braking determined via the period 
derivative P are well measured and constrained by observations, the topology and strength of the magnetic field anchored 
in the neutron star is much less well known, with large uncertainties. The magnetic field can only be inferred by some 
assumptions like for instance about the spin-down luminosity induced by magnetodipole braking. Magnetic field estimates 


from the cyclotron line emission is a notable exception where severe constraints can be put to good accuracy aTruemper et al 
19781') . 


With no a priori knowledge about the magnetic field configuration, the most simple structure is assumed namely a dipole 
field. A rich literature exists on the consequences of this field around neutron stars. The dipole field is almost exclusively used 
to describe the magnetosphere of pulsars and magnetars. Although such field can sometimes fairly explain light-curves from 
the pulsed radio and high-energy emission, it sometimes completely fails to account for the polarization properties of other 
radio pulsars. It is hard to believe that multipoles are not anchored in the neutron star crust as such fields are already present 
in main sequ e nce s tars such as the Sun. 

Deutschl (I1955T I was the first to compute the electromagnetic field around a rotating dipole, taking into account the finite 
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size of the star. His solution is often quoted to explain the magnetic dipole radiation losses from neutron stars. Neverthe¬ 
less, higher multipoles can also contribute to the spin-down luminosity, to a fraction depending on the ratio R/r l and on 
Bmuitipoie/Bdipoie where R is the neutron star radius and n. = c/fi the light-cylinder radius with c the speed of light and 
= 2n/P the star rotation rate. Effects of higher multipole moments can be drastic especially close to the neutron star 
surface where they could be dominant. 

However accounting for neutron star braking only from magnetodipole losses is not possible due to an observed braking 
index n always less than 3 which contradicts the point magnetodipole losses. Some part of the Poynting flux has to be 
explained by another multipole or for instance by a particle outflow like a relativistic wind of electron/positron pairs, leading 
to a braking index of n = 1 . Let us briefly mention some previous works about multipolar magnetic Helds in neutron stars. 
Blandford fc Romani ( 198gli discussed the consequences of a non-dipolar magnetic field on the spin-down evolution of pulsars. 


They took in to account a non-c onstant in time proportionality factor for the deceleration with application to pulsars available 
at that time. I Chen fo Lil (]200£) explained the braking indexes less than three by an evolving magnetic dipolar moment or the 
presence of a fall-back disk. (Yue et alj (120071 ') investigated a two component model between the pure vacuum rotating dipole 
and the pure polar cap relativistic outflow assuming a constant potential drop. They found a braking index always between 
the two extremes given by n = 1 and n = 3. From an observational point of view, radio emission of pulsars usually requires the 
formation of vacuu m gaps. It seems that for some pulsars, the surf ace magnetic field needs significant multipolar components 
( Gil fe Mitral 2001 1. This is especially true close to the polar caps ( Gil et al. 2002a 3). 

Investigations of t he effect of mag netic multipoles on the neutron star shape and its surrounding have also been done 
during the last decades. Roberts ( 1979ll computed the electromagnetic multipolar fields to deduce the electromagnetic recoil of 
a freshly borned neutro n star. He also computed the off-centred magnetic multipoles and brought some higher order corrections 
to the Deutsch s olution. Krolik (1 1991 1 showed that for millisecond pulsars constrains can be put on the strength of the magnetic 
multipole fields. Mitra et al.l ( 1996 1 studied the decay of multipolar ma gnetic components in neutron stars and did not find any 
significant evolution able to modify the radio pulse profile. Recently, Mastrano et al. ( 2013l l looked at the dipole-quadrupole- 
hexapole magnetic field induced deformation of a neutron star opening up the possibility for emission of gravitationa l wave s 
showing various patterns and discussed implications for the braking index of magnetars. Asseo fe Khechinashvili ( 2002l l 
investigated the consequences of the presence of multipolar fields in pulsar magnetospheres close to the surface. Special 
attention was drawn on pair cascades with application to the sub-beams of PSR B0943+10 and some sub-pulse drifting 
phenomena. However, the y mainly focused on ax i symm etric configurations that do not permit pulsed emission. Pair cascade 
was also the main topic of Harding fe Muslimovl (2011) who found that a lready in the disto r ted d ipole field the accelerating 
potential increases significantly. A simi l ar inv estigation was performed by Kantor fe Tsvgan ( 2003h with application to crust 
heating and electron extraction. Arond ( 1993l l looked at the spin up of neutron stars due to accretion and found a constraint 
on the multipole components in the magnetosphere, claiming that it will always be dominated by the dipole. Recently, 


Bonazzola et alj (1201.51 1 computed the solutions to the multipole field in vacuum using two scalar fields that automatically 


satisfy the divergencelessness constrain for the magnetic field. Our approach with vector spherical harmonics is similar to their 
scalar fields although that our formalism is more handy to solve for the electromagnetic field components in vacuum. They also 
briefly discussed the consequences on pulsar electrodynamics and pulse shape. Let us also mention that general-relativistic 
mu ltipolar electromagnetic fiel d decomposition around a rotating neutron star has been carried out more than a decade ago 
bv IRezzolla fe Ahmedovi (120041 '). They estimated the corresponding spin-down luminosity for a rotating magnetic dipole and 
showed it to be 2 to 6 times larger than the often quoted Newtonian limit. 


Alvarez fe Carraminana ( 19981 ') studied the evolution of pulsars according to the spin-down torque produced by a monopole 


(wind outflow), adding a dipole and a quadrupole co ntribution. Lipunov ( 1983h even discussed the possibility to detect 
multipole radiation from neutron stars. Shi fe Li ( 19981 ') studied the gamma-ray emission coming from the eight gamma-ray 
pulsars known at that time. Curiously they found that the quadrupole fields was as strong as the dipole field suggesting a 
relation between high-energy emission and magnetic topology. 

Understanding the very many classes of neutron stars probably requires a deep investigation of the exact topology of 
the multipolar fields to unveil the physics of neutron star magnetospheres. The clue resides in the relative strengths of each 
multipole moment. In this paper, we study in details the properties of the first four multipole moments, namely the magneto¬ 
dipole l = 1, the magneto-quadrupole l = 2, the magneto-hexapole l — 3 and the magneto-octopole l = 4 losses. In 2] we give 
exact close analytical expressions for the electromagnetic field of a rotating multipole and for the corresponding Poynting flux 
in terms of spherical Hankel functions. We then specialize our treatment to the dipole, quadrupole, hexapole and octopole 
in Sj3] In the appendix [B] we also give simple approximations for those multipoles in the limit of a vanishing radius of the 
star, the so called point multipole approximation, valid to good accuracy for all stars with R ^ n, which is the case even for 
millisecond pulsars (for which R/ri, < 0.1). The consequences of such fields are then investigated. The associated Poynting 
flux and braking indexes are given in [21 Upper limits for the magnetic multipole components are given according to the 
magneto-multipole losses formula. The polar cap geometry is strongly influenced by the non-dipolar components as shown 
in m The pulsed high-energy emission within the magnetosphere, along the last closed field line (slot gap model) shows 
drastically different properties compared to the dipole as discussed in m Conclusions and ongoing work are drawn in an 
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Multipolar electromagnetic fields 3 


2 EXACT MULTIPOLAR ELECTROMAGNETIC FIELDS 


Finding the multipolar electromagnetic field of a rotating magnetized sphere can be reduced to a set of homogeneous wave 
equations for scal ar quantities associated to the electric and magnetic field. The general fr amewo rk in general relativity 
has been given in Petri ( 20131 '). Here we specialize the discussion to Minkowski space-time ( Petri 20121 '!. avoiding frame¬ 
dragging effects which would generate higher multipoles from lower ones. All the details about the derivation of th e partial 
differential equations recalled below and the associated expansion method are taken from chapter 9 of ljacksonl (120011 ) dealing 
with radiating systems and multipolar fields. However, our definitions do not strictly follow the usual conventions taken for 
quantum mechanics. We give exact analytical solutions for any multipole moment, specifying boundary conditions on the star 
being treated as a perfect conductor. 


2.1 General solution 


Let us assume that the neutron star interior is filled with a constant magnetic field dragged into solid uniform rotation by the 
star. The magnetic field can then be cast into an expansion according to the vector spherical harmonics such that in spherical 
coordinates ( r , i?, <p) 

i 


Bin(r,0,tp,t) = ^ X! ( V x 


(1) 


The exponential term emphasizes the solid body rotation of the magnetic field structure as seen by a distant observer 

in its inertial frame. We recall that the azimuthal term e lrnv can be factored out of such that each mode m behaves as 
i) corresponding indeed to corotation at the neutron star speed fi. For r ^ R, the functions /; B m (r) are prescribed 
by the particular magnetic field configuration inside the star. These are known functions and will impose the boundary 
conditions on the neutron star surface. Note that the dependence on time is entirely contained in the harmonic terms, the 
complex exponentials it is u nderstood th at the physical solution corresponds to the real part of B ln . The vector 

spherical harmonics <I>(, m are defined in Petril ( 20131 ). To the lowest order in the expansion, there is no monopole field therefore 
no l = 0 term. We are interested in the exact multipolar electromagnetic field induced by a rotating neutron star in a flat 
vacuum space-time. The divergencelessness constraint on the electric field and magnetic field can be conveniently expressed 
in terms of vector spherical harmonics Indeed the electromagnetic field outside the star is also developed on to vector 

spherical harmonics such that 

i 


D ou t()",#,<£,£) = ']T ^2 (v x [/i D m (r,f)§i, m ] + 


(2a) 


Bout(r,tf,v3,£) 5Z ( v x [/i B m(r,£)$!,m] + Si B m(r,f)$i, m ) . (2b) 

1=1 m= — l 

We are looking for stationary solutions expressed in the observer frame. This is achieved by a temporal dependence of 
and gi,m given by Thus the functions gf > m and gf rn are related to the function /; B m and / B m according to Maxwell 

equations by a simple linear scaling such that 


D 

9l,rr. 

B 

9l,rr 


= + i Eo mft „ 
= - i go m fi fi 


(2c) 

(2d) 


We recall that in flat space-time vacuum the constitutive relations are D = eqE and B = go H. For a general multipolar 
electromagnetic field, the functions {/ B m ,/i D m } are separately solutions of an Helmholtz equation given by 


; dr(rfi, m ) - 


1(1 + 1) 


fl,m + km fl,m = 0 


where we introduced the wave-number k m = mk and 


k_ — _ — 

C TL ' 


( 3 ) 


( 4 ) 


In such a way, we decoupled the electric part /; D m from the magnetic part fi m . This decoupling arises from the definition of 


the vector spherical harmonics «1?( 


being eigenfunctions of the angular Laplacian differential operator in the sense that 
1 d . „ d 1 d 2 " 


d_ . d_ _ 

sin $ dP Sln <9$ sin 2 P dip 2 


&l,m = -1(1+ 1) $l,r, 


( 5 ) 


Equation ([3]) and its general solution can be found in any textbook such as for instance in I Jackson! (120011) . The solution 
corresponding to a spherical outgoing wave, that with m > 0, is 

( 6 ) 


fern — ^7 (^m r ) 
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where h[^ ( x ) is the spherical Hankel function of order l as defined by Arfken fc Weber ( 2005l l . The constant ai :Tn is determined 
by the boundary conditions on the star, namely continuity of the radial component of B and continuity of the tangential 
component of D. We will come to this point later again. Moreover the axisymmetric modes m = 0 do not propagate. They 
decay to zero at infinity according to a power law in radius given by 

/i,o = «i,or" (1+1) . (7) 


ai : o are also constant to be determined by the boundary conditions that, on the neutron star surface, enforce continuity of 
the radial component of the magnetic, implying continuity of the functions / ( B m at the stellar surface. For the electric field, 
we must satisfy for any m 


dr(r / ; D m ) =£ 0 m [V(J+1)(Z-1) Jl,m ft 1, B - Vl (l + 2) Jl+l,m + 


( 8 ) 


where Jz, m = see Petri ( 2013l l. These relations fully determine the electromagnetic field outside the star for any 

multipolar component. Indeed, the general solution reads 


OO _ OO l 

Pout (r, if, t) = y V x [q° 0 -4^] + y (v x [af }Tn h\ 1} (km r) $z, m ] + ie 0 m H a B m h, (1) (fc m r) e~ 


1=1 m= — l,m^ 0 


OO ^ OO L 

B out {r,&,ip,t)= y^V x [a B 0 + y y (^V x [a® m /i, (1) (A m r) $i, m ] - i m fl a° m ftj 1 ’ (fc m r) j e~ 


1=1 m= — I,m^O 


(9a) 


(9b) 


where {aP m , a B m } are constants depending on the boundary conditions imposed on the stellar surface. Note that the presence 
of a B m terms in the magnetic field and of af m terms in the electric field does not contradict the fact that equations @ are 
decoupled. The coupling is due to the gi iTn functions as defined in equations (I2cl) . (l2dl ). Any particular solution to the field 
equation requires the determination of these constants {aP m ,af m }. First the continuity of the radial component of B for 
m > 0 gives 


_ fl,m(P) 
0*1,111 7TC~ “ 

hY\k m R ) 

and for the axisymmetric cases m = 0 

alo = R l+1 ft(R)- 


( 10 ) 


( 11 ) 


Second, the continuity of the tangential component of the electric field gives for the axisymmetric and non axisymmetric 
coefficients respectively, 

UP 0 = -eo VR l+2 [y(I + 1) (l ~ 1 ) Ao f?-i,o(R) -#T2)li +1 ,o/| B +1 ,o(fi)] (12a) 

af,m d r (rh[ 1) (k m r))\^ ii = £oRn [V(Z + 1) (Z - 1) Ji, m fti,m(R) ~ y/l (1 + 2) Ji +1 , m f?+i,m,(R)\ (12b) 

where derivatives are evaluated at the stellar surface r = R. This represents the most general expression to the vacuum 
multipolar electromagnetic field outside the neutron star. Knowing the radial component of the magnetic field on its surface, 
we are able to compute straightforwardly the full electromagnetic field in vacuum. In what follows, we specialize these solutions 
to some practical situations including low order multipoles like the magnetic dipole, quadrupole, hexapole and octopole. We 
also give useful and handy compact expressions in the limit of vanishing radius of the neutron star: we refer to it as retarded 
point multipole field, see 331 and the appendix [B] Before, we show how to compute the exact multipole Poynting flux for any 
mode ( l,m ). To this end, we look for the asymptotic expression of the electromagnetic field in the wave zone. 


2.2 Wave zone and Poynting flux 


In the wave zone, the above expressions can be drastically redu ced by the fact that the spherical Hankel functions behave 
asymptotically as h^ 1 (x) « (— i) l+1 e lx /x ( Arfken fe Weberll2005l l. Neglecting the axisymmetric mode decreasing much faster, 
like r - ^ i+i P the electromagnetic field becomes in the limit r r l 

i (fc m r—m f2 £) 


Bw= y -(-*) 
D w= y -(-*) 

Z^ 1,771^0 

=£o cB w A n. 


i e 


i e 


km r 

i (fc m r — m f2 t ) 


km T 


CLl,m ^Z,77i 0 ^Yl CLi^m ^Z,??!^ 

{^km &l,m ^ Z ,771 £0 Tri ^ ^Z,77i ^Z,77i^ 


(13a 


(13b) 
(13c) 
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l/m 

1 

2 

3 

4 

Normalization 

1 

2 

3 

1 

32 

to 

64 

1&4 

729 


a 4 

a 6 

a 8 

4 

m 

99225 

4§S§6 

99225 

13^-6 

1225 

1048576 

11025 

a 10 


Table 1. First order expansion of the function Si m in terms of a — R/r l- The first column indicates the multipole order l and the 
first line indicates the azimuthal mode m with the restriction m ^ l. The normalization is proportional to a 2 * - * -2 oc fi 2i + 2 , typical for a 
/-multipole. 


Equation (I13cl) shows that the solution behaves as a monochromatic plane wave propagating in the radial direction n = e r at 
frequency Q. It is derived from the identities 


n A ^l } m = 

n A 


The time averaged Poynting flux is therefore 


S = 


Dw A b; 


(14a) 
(14b) 

2,0,0 (15) 
where BJ, is the complex conjugate of B w . Integrating the radial component of the Poynting vector along the solid angle we 
get the power radiated, using the orthonormality of the vector spherical harmonics, such that 

L= f S T r 2 dfl = —— (\a?,rn\ 2 + UoC 2 \afm.\ 2 ) ■ (16) 

■* n ^ 0 1^1,m^O V 

The spin down luminosity L is independent of the radius as it should from the energy conservation law. Equation m 
represents th e most general expression for the magneto-multipole losses from an arbitrary multipole magnetic field, see for 
instance also Jackson! ( 200 il l. 


2.3 Exact solution for one multipole 


It is instructive to find the exact solution for a particular multipole held with fixed numbers ( l , m). Let us assume that inside 
the star, the magnetic held is solely represented by the function /j B m (r). What should then the electromagnetic held be outside 
the star? The only non vanishing magnetic held coefficient is given by equation Col) or equation CD if m = 0. Moreover the 
two non-vanishing electric held coefficients are (if l = 1 only one solution exists, see the dipole case below) for m > 0 


aj+i.m d r (rh[f! 1 {km.r))\ = e 0 RQ. ^l (l + 2) ffm(R) (17a) 

I r=R 

dr(r hfffkm r))| = -e 0 R Q (l - 1) (l + 1) Ji, m ff m (R) (17b) 

I r=R 

and for the axisymmetric case we hnd 

0 + 1) aP+i,o = -eo R l+3 V Vl (l + 2) Ji+i,o ffo(R) (18a) 

(l ~ 1) aP_ li0 = eo R l+1 SI y/(l- 1) (Z + lj Ji, o ffo(R)- (18b) 


We conclude that 
associated to this 


the solution is fully specihed by the three constants of integration (a B m , af + l m , ap_i,m)- The Poynting hux 
particular solution is for m > 0 (it vanishes for m = 0) 


Pl.m = 


B i2 
^ l.m | 


Sl,m — 


2/i0 

c\ffm(R)\ 

2 ,o 

1 


+ Po C 2 (| a B-l,m | 2 + | a Pn,m| 2 )j 


\h { f\kmR)V 


w_ 

rt 


lil + 2) 


+ 


(i) 


\drirh\fpkm r))\ 2 R \d r (r h\\(k m r))|| 


(19a) 

(19b) 

(19c) 


For the special case l = 1, the constant ap_i,m does not exist. For later purposes, it will be useful to get the series expansions 
to hrst order of in terms of a = R/r l- For completeness, we give them in Table C] Actually they represent within a 
multiplicative factor the Poynting hux for the retarded point multipole held. For concreteness, we switch now to explicit 
application of low order multipole solutions in order to get more physical insight into their properties. For the remainder of 
this paper, we focus on some illuminating cases such as the low order multipoles: dipole, quadrupole, hexapole and octopole. 
To end this general discussion about the computation of multipolar helds, we compare our results with the recent work done 
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by Bonazzola et al. 1 20151b It is shown in the appendix [X] that we get exactly the same expressions as Bonazzola et al. 1 2015h 
except for an error in the sign of their electric field for the axisymmetric multipoles. 


3 EXPLICIT LOW MULTIPOLE SOLUTIONS 

Multipolar magnetic fields are important not only for neutron stars but also for main sequence stars, or any other star 
possessing a non negligible magnetic field. To study quantitatively the influence of those components, it is illuminating to 
get simple closed analytical expressions for those fields. Thus in this section, we give the full exact solutions for the dipole, 
quadrupole, hexapole and octopole fields. For practical applications, we also give the approximation for a retarded point 
multipole useful and very accurate for slowly rotating stars as those we study here: by slowly rotating we mean R r l. 


3.1 The magnetic dipole l = 1 


We start our discussion with the well known magnetic dipole. The ease and power of our description is demonstrated by 
application to the dipole field. Introducing the vector spherical harmonics expansion and two free parameters {qi, 0 , 91 , 1 } 
related to the strength of each azimuthal mode m, the quasi-static rotating dipole frozen into the neutron star is conveniently 
written as 


B= BR* Re Vx 


$1.0 $1,1 -iOf 

<7i,o - 3 —h < 71,1 —— e 


r , - . r , , (20) 

The symbol Re means the real part and the two constant coefficients { 51 , 0 , 91 , 1 } are traditionally expressed in terms of the 
magnetic obliquity x, angle between magnetic dipole vector and rotation axis, such that 

/ 8 7T 


<71,0 = — 


<7i,i 


cosx 


16 7T 


smx- 


(21a) 


(21b) 


Thus x = 0° gives the axisymmetric mode m = 0 alone whereas x = 90° give the m = 1 mode alone. Expression flU} is 
only valid close to the star where retardation effects are negligible. It corresponds to the quasi-static limit, the lowest order 
approximation in R/r l. Actually it also represents the boundary value imposed by the rotating neutron star. Solving the 
Helmholtz equation m for the unknown expansion functions {/i D m ,/j B m }, the rotating magnetic dipole inside the star is 
advantageously expressed by only two non-vanishing functions given by 


fro = ~\l — 


8^t B R 3 cos x 


a _ ,16 7 1 B R 5 sinx 

■' 1 ’ 1 “ V ~~3 r 5 ' 

The boundary conditions impose continuity of {/go, / B 1 } at the stellar surface and for the electric field we get 

9r(r / 2 ,o) = -y= £0 RVL /go 

f Eoflfi^i. 

5 


d r {rf 2 ,i) = 

The full solution in vacuum outside the star therefore becomes 


B 

a i,o 


B 

a l,l 


-\J^BR 3 cosx 


D 

a 2,0 ^ 


D 

a 2,l — 


16 7t B R sin x 

3 h ( i\kR) 

/ 8 7T 5 

—— £0 12 B R cos x 
15 


16 7t ^ VLB R 2 sin x 
d r (r h^\k r))\ 


■£o 


(22a) 

(22b) 

(23a) 

(23b) 

(24a) 

(24b) 

(24c) 

(24d) 


which is nothing else than Deutsch solution ( Deutsch 19551) . It is easy to check that this configuration with only one free 
parameter corresponding to the inclination angle of the magnetic moment x leaves the total magnetic energy outside the star 
constant. Indeed the total magnetic energy amounts to 


Wmag = 


mi 


dV = 


B 2 R 2 
4 Mo 


(2 < 71,0 + <7i,i) = 4 7t 


B 2 R i 
3 Mo 


(25) 


© 2013 RAS, MNRAS OOO.IT11291 



































Multipolar electromagnetic fields 7 


independently of x- dV stands for the integration in the whole three dimensional space outside the star. The total magnetic 
energy in vacuum outside the star will be used to scale the different multipole moments with respect to the dipole one. 


3.2 The magnetic quadrupole l = 2 


Let us perform the same expansion to the magnetic quadrupole such that it can be expressed inside the star by 


B = BIT Re Vx 


3 ? 2,0 3 ? 2,1 -in t , ^2,2 -2ifit 

92,0 —,-b 92,1 —e + 92,2 —e 


( 26 ) 


The boundary conditions impose continuity of {/® 0 , f?,i, f 2 , 2} a t the stellar surface and for the tangential component of the 
electric field we get 


9 r(r /i,o) — £0 RQ f 2 ,0 


drif / 11 ) = ~\ -E 0 Rttf 2 ,i 
V 5 

dr(rf^ 0 ) = 6^eoR^ff, 0 
d r (r so Rnf? A 


dr(r fi'T) — 2 y ^ £ o RQ f 2 , 2 - 

Solving again Helmholtz equation ©, the solution in vacuum is represented by 


D T-) T-) 4 

“2,0 — 92,0 B H 


a 2 ,l — 92,1 


a 2,2 = 92,2 ■ 


BR 


R) 

BR 

h ( fi(k 2 R) 


n?o = — 7= £0 SIB R 4 92,0 
v 5 


“1,1 = -\ M £ o ■ 
5 


B R“ 92,1 


d r (r h^fiki r))| 


03,0 = — 2 \j — eqQ.BR 92,0 


D 

“3,1 — 


fi 35 d r (r r))| 


fl B R 92,1 


H B R 2 Q 2,2 


“3,2 ^ 2 \ j - £ 0 ' 

7 d r (rh 3 (fer)) 


( 27 a) 

( 27 b) 

( 27 c) 

( 27 d) 

( 27 e) 

( 28 a) 

( 28 b) 

( 28 c) 

( 28 d) 

( 28 e) 

( 28 f) 

( 28 g) 

( 28 h) 


We have to fix the three constants {92,0,92,1,92,2} in a judicious way to compare results with the dipole. To do this, we 
compute the total magnetic energy outside the star from its quadrupole moment and find 


mi 


_ B R tr, 2 .2 | 2 \ 

dV — — -( 2 q 2 0 + 92 1 + 92,2)- 

2/ro 


( 29 ) 


Assuming a constant magnetic energy, irrespective of its geometry, the coefficients {92,0,92,1,92,2} can be related by 


92,0 = 

92.1 = 

92.2 = 


cosxi 


n . 

sin xi cosx '2 


sin Xi sm x2- 


( 30 a) 

( 30 b) 

( 30 c) 


where Xi £ [ 0 , 7t] and X2 £ [ 0 , 2 7t] are two angles specifying the particular geometry of the quadrupole magnetic field. With 
this normalization of {92,0,92,1,92,2}, the magnetic energy is the same as for the dipolar field. The different modes are easily 
separated by taking xi = 0 ° for m = 0 whatever X2, (xi, X2) = ( 90 °, 0 °) for m = 1 and (xi, X2) = ( 90 °, 90 °) for m = 2 . 
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3.3 The magnetic hexapole l = 3 

We show how our approach is straightforwardly extended to hexapole or even higher multipolar fields. Indeed we 
the magnetic hexapole existing inside the neutron star such that 


B = BR b Re Vx 


93,0 ■ 


3?3,0 


+ 93,1 


‘JWi -i 


+ 93,2 ' 


^3,2 -2iQt 


+ 93,3 ' 


-3,3 -3int 


The boundary conditions enforces continuity of {/ 30 , / 3 ,i> / 3 , 2 i ^ 3 , 3 j- at the stellar surface and 


dr(rf^ 0 ) = -6 ] j^eoRnf^ 0 

dr(r /°i) = — —^== so RQ ff }1 


d r {rf 2 , 2 ) = -2 x /-£ 0 RQf^ 2 


dr{rf 4 ,o) =4:\J — £ 0 R fi/ 3,0 

d r (rf° r) = -^=E 0 Rttfi A 


dr(r f?, 2 ) = 2 x /- eo RV tf 2 


dr(r / 4 , 3 ) -J-e 0 RQ, f 33 


Computing the exact solutions to the Helmholtz equations © leads to 


« 3 ,o = 93 ,o B R° 
b BR 

a 3,l = 93,1 


h£\ki R) 
b BR 

a 3,2 = 93,2 -Try-- 

h 3 (k 2 R) 

B BR 

a 3,3 — 93,3 ' 


h^(k 3 R) 


« 2 ,o = 3-\/ — 53,0 


D 

a 2,l — ~ 


V35 d r ( r h£\kir))\ 


Q B R 2 53,1 


fl B R 2 53,2 


a 2,2 — -2 \l - £0 ■ 

^ d r (rh 2 \k 2 r))\ 


D 

a 4,0 


D _ 5 

G .4 I - - — £0 


— eoQB R 7 53,0 
fl B R 2 53,1 


d r (r h^\kir))\ 


LJ O / 

a 4 o — ^ \ — £0 


Q B R 2 q3,2 


d r {rh\ \k 2 r)) 


D 

a 4,3 


5 fl B R 53,3 

^ dr{rh ( i\k 3 r))\ 


To compare with the dipole field, we compute the total magnetic energy as 

3 B 2 R 3 


Wmag = 


///£ 


dV = 


4/io 


(2 53,0 + 93,1 + 93,2 + 93 , 3 ) 


start with 

(31) 

(32a) 

(32b) 

(32c) 

(32d) 

(32e) 

(32f) 

(32g) 

(33a) 

(33b) 

(33c) 

(33d) 

(33e) 

(33f) 

(33g) 

(33h) 

(331) 

(33j) 

(33k) 

(34) 
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Assuming a constant magnetic energy in the field, irrespective of the geometry, the coefficients are related by 



2 V2tt 

93,0 = 

3 cos Xi 


4 . 

93,1 = 

3 sin xi cos X2 


4 \fri . 

93,2 = 

—— sin xi sin X2 cosx3 


4 y/n . 

93,3 = 

—— smxi sin X2 smx3 


(35a) 

(35b) 

(35c) 

(35d) 


where xi, X 2 £ [0,7t] and X 3 £ [0, 2 7t] are three angles specifying the particular geometry of the hexapole magnetic field. With 
this particular normalization, the magnetic energy is the same as for the dipolar field. The different modes are easily separated 
by taking xi = 0° for m = 0 whatever (X 2 ,X 3 ), (XI 1 X 2 ) = (90°, 0°) for m = 1 whatever X 3 , (xi,X 2 ,X 3 ) = (90°,90°,0°) for 
m = 2 and (x'1jX 2,X3) = (90°, 90°, 90°) for m — 3. 


3.4 The magnetic octopole l = 4 

To finish this discussion about the low order multipoles, we give the exact solution to the octopole field. Inside the neutron 
star, the magnetic field is described a priori by a general octopolar expansion such that 

| (36) 

The boundary conditions for the non vanishing electric field components are listed below, 


B=BR 6 Re(Vx 


4 *4,0 3?4,1 -int . 4>4,2 _ 2 int . 4*4,3 . 4>4,4 

?4,o —;-h 94,1 —r— e + 1 / 4,2 —r— e + 94,3 —r— e + 94,4 —r— e 


dr(rf£ 0 ) = -4 ] j^£oRnff y0 

(37a) 

dr(r /°i) = —-^= eo RH f 4,1 

(37b) 

dr{rfl 2 )=-2^e 0 RBlfl 2 

(37c) 

dr{rfZ 3 ) = -^soRSlfZ 3 

(37d) 

dr(vfZo) = W ] j^£ 0 RnfZ 0 

(37e) 

dr{rfZ 1 ) = -j=e 0 RnfZ 1 

(37f) 

dr{rfF A )=2^e 0 R^fZ 2 

(37g) 

dr(rf?, 3 ) = 8^e 0 RnfZ 3 

(37h) 

dr(rfF A )=2^soRnf? A - 

(37i) 
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Solution to the Helmholtz equations Q therefore leads to the complete solution in the form of its constant of integration 
given by 


“ 4,0 — 54,( 


a 4,l — 54,1 ttt 


BR 


°4,2 — 54,2 —yjy 


PPik! R) 

BR 


hi l> (ka R) 


a 4,3 — 54,3 —(jy 


BR 


h^(k 3 R) 


a 4,4 — 54,4 —yjy 


BR 


h { P{k 4 R) 


D 

a 3,0 


— E 0 nBR 6 q 4:0 

Q B R 2 ^ 4,1 

- —J= fco-—-- r 

V 7 d r (r /13 (fci r)) 


a 3,2 


= -2 


£0 ■ 


n b r? ^ 4,2 


“3,3 — o £0 


7 dr(r h^\k 2 r)) 

H H R? 54,3 
^ 9r(r - ft3 1) (fe3r))| 

I r=f 

= — 2 7 ^ £qQ. B R s 54,0 

B jR 2 $ 4,1 


r=R 


D 

a 5 ,l 


D 0 / J- 4 * 

a 5,2 — ^ A/ TT £0 


^ r (r r))| 
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Q B R 2 Q4,2 


-CO -j- 

dr{r hl 1 \k 2 r))\ 


r=R 


a 5 ’ 3 - 8 V « £ ° 


H 1/4,3 


^3 d r {r hf$\k 3 r))\ 


\r=R 


a 5,4 


D =2j^-£ 0 


H B R“ 54,4 


The total magnetic energy outside the star is is 


H ^r(r /l 5 ^(/C 4 r))| 


r=H 


—!!!&"-** 


2 2 2 2 2 
(2 54,0 + 54,1 + 54,2 + 54,3 + 54 , 4 ) 


"mag - / / / 7 ;- av — - IA 54,0 T 54,1 “T 54,2 T 54,3 ~r 54,41 

J J J 2 Ho fio 

Assuming a constant magnetic energy in the field, irrespective of the geometry, the coefficients are relate 


(38a) 

(38b) 

(38c) 

(38d) 

(38e) 

(38f) 

(38g) 

(38h) 

(38i) 

(38j) 

(38k) 

(381) 

(38m) 

(38n) 

(39) 


54,0 = 


54,1 


■ — 


54,3 = 


54,4 = 


cosxi 


sill Xl COS X2 


sm xi smx 2 cos X3 


sm xi smx 2 smx3 cosx'4 


471 . 

— sm xi sm X 2 sm X3 sm X4 


3 

471 

3 

471 

3 

471 


(40a) 

(40b) 

(40c) 

(40d) 

(40e) 


where xi, X2, X 3 £ [0,7i] and X4 £ [0, 2 71] are four angles specifying the particular geometry of the octopole magnetic field. With 
this particular normalization, the magnetic energy is the same as for the dipolar field. The different modes are easily separated 
by taking xi = 0° for m = 0 whatever (x 2 , X3, X4), (x i,Xa) = (90°, 0°) for m = 1 whatever X 3 , X 4 , (xi, X2, X3) = (90°, 90°, 0°) 
for m = 2 whatever X 4 , (xi, X 2 , X3, X4) = (90°, 90°, 90°, 0°) for m = 3 and (xi, X 2 , X 3 , X 4 ) = (90°, 90°, 90°, 90°) for m = 4. 

More handy expressions for these multipoles are given in the appendix m where approximate solutions in the limit 
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of a point multipole are shown (in the limit R —> 0). The retarded point multipoles represent very accurate and simple 
analytical expressions to study the effects of non dipolar fields on the structure of pulsar magnetospheres (or other non 
compact objects), its polar caps and its emission geometry. The detailed consequences of such fields, for instance on the radio 
polarization properties, requires a thorough analysis which we leave for future work. However, to get a flavour of multipolar 
effects, for the remainder of the paper, we will show the influence on the Poynting flux and related braking index as well as 
on the polar cap geometry and high energy emission phase diagrams. 


4 POYNTING FLUX AND BRAKING INDEX 

As a diagnostic of the consequences of the presence of multipole fields, we compute the exact Poynting flux of each multipole 
labelled by the mode (l, m) taking into account the finite size of the star. The braking index n is an interesting related 
quantity which describes the efficiency of electromagnetic radiation by the law summarized in Cl = — K QL where A is a 
constant depending on boun dary conditions on the neutron star. For magnetic multipolar point sources of order l, we know 
that n = 2 1 + 1 (IKrolikfll99lh but in general it can differ from this fiducial value if the size of the star is taken into account. 
Formally, the braking index follows from the luminosity L (or the power radiated) by the relation 

9 dL 


n = -1 

L dn 


which is derived from the definition of the braking index as 


and from the spin-down luminosity given by 


no. 

7F 


l = -i fin 


(41) 


(42) 


(43) 


where I represents the moment of inertia of the neutron star. Note the minus sign in order to have L > 0 because braking 
implies 9 < 0. Expressed in terms of the dimensionless parameter a = R/r l we get 


_ a dL dlnL 

L da d In a 


(44) 


This is the general formula to compute the braking index in any case, knowing the luminosity of the star with respect to the 
spin normalized by the parameter a. Useful approximate expressions for the radiated power L and the braking index n are 
given in the following paragraphs starting from the dipole field and up to the octopole. 


4.1 Dipole 

First, we recall the exact Poynting flux obtained from Deutsch solution and compare it to the point dipole approximation. 
The magnetodipole point source losses are 


j- point _ 8 71 

dip " W 


9 4 B 2 R 6 


sin 2 x 


which can differ significantly from the true dipole losses generated by a star of finite size, according to the formula 

Ldip _ 1 _3a 4 _ 


t point 
-^dip 


(45) 


(46a) 


i 2 + 1 5 (a 6 - 3a 4 + 36) 

rs 1 — a 2 + o(a 2 ) (46b) 

Equation (I46al) is consistent with the expression given by Melatosl (1997). From our general treatment of the multipole fields, 
we would write it as 

1 3 K* 1 


T 8 7T r>2 j-,2 • 2 

Ldip = -—cB R sin X 
3/r o 


_ 3 Rfi _ 

\hfi(k R)\ 2 5 \d r (r fi'fi (k r))| 2 


flj 


(47) 


which is exactly the same as equation (I46al) . The spin-down luminosity is always less than the point dipole source, its variation 
with respect to the stellar rotation speed is shown in Fig. [TJ assuming that only the mode m = 1 is present (x = 90°). Because 
for all known pulsars, even for millisecond ones, we have a < 0.1, the approximation to second order in a, as given by 
equation (I46bl) . is always sufficient to compute accurately the spin-down rate. The luminosity decreases monotonically with a 
reaching only 52% of the point dipole source for a = 1 although this rotation rate becomes unrealistic. Nevertheless it shows 
that the spin-down luminosity slightly changes with respect to 9 and thus, rigorously, the braking index will deviate from the 
often quote value of three. Indeed, this braking index derived from equation (1441) is given by 

2 6 (a 4 - 45) 6 (a 4 - 36) 


n = 1 + 


J 2 + 1 2a 6 - 3a 4 + 45 

3 — 2 a 2 + o(a 2 ) 


a 6 - 3a 4 + 36 


(48a) 

(48b) 
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R/r L 
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Figure 1. Spin-down luminosity of the rotating magnetic multipole according to the ratio R/r l. The contributions are separated into 
single components (Z,m). The dipole is shown in red, the quadrupole in green, the hexapole in blue and the octopole in magenta. Note 
that each luminosity is normalized to its value for R = 0. 


The braking indexes given in this paper have been com puted wit h the formal calculator Mathematica. It is straightforward 
to check that eq. (I48al) is exactly identical to eq. (B4) of Melatosj (1997) in the special case of p = 1. Indeed, in our picture, 
the empty space corresponds to the region outside the star therefore using his parameter defined in his equation (6) we have 
x v = RQ/c thus x v /x v = Cl/n implying p = 1 according to his definition. The braking index strictly equals three only for 
the retarded point dipole. It is always less than three although not significantly except for fast rotating star for which a < 1. 
For instance, for millisecond pulsars with a = 0.1, the braking index is still as high as n = 2.98. Its variations are shown in 
Fig. [2] Starting from n = 3 for a = 0 it decreases to n = 2.18 for a = 1. 

We must conclude that the Poynting flux and its braking index are subje ct to var i ations depending on the location of 
the inner boundary with respect to the light-cylinder as already mentioned by Melatosl (1997). This tells us that a braking 
index less than three can be explained by artificially increasing the size of the neutron star to such a point that it reaches 
almost its light-cylinder. This could happen for instance if the magnetosphere is only partially filled with electron/positron 
pairs, mostly in corotation with the star in a kind of dome+torus shape but also with an outer disk in ov errot ation. Such 
plasm a equilibria are known for a long time and referred as electrospheres by Krauge-Polstorff fc Michel ( 19851 ) and Petri et~al1 
( 20021 ). This argument is similar to the corotating inner magnetosphere of Melatosl " 1997 1 expect that in the latter model, 
the boundary is kept spherical whereas in the former, the boundary surface is much more complicated. The switch between 
a vacuum magnetosphere and an electrosphere leads to a switch in the P oynting flux and therefore could explain the switch- 
on/switch-off modes of several intermittent pulsars ( Kramer et al. 20061 ) although a quantitative analysis would require an 


accurate value of the spin-down power of the electrosphere, which unfortunately is not yet accessible. 
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Figure 2. Braking index of the rotating magnetic multipole according to the ratio R/r l. The contributions are separated into single 
components (l, m). The dipole is shown in red, the quadrupole in green, the hexapole in blue and the octopole in magenta. 


4.2 Quadrupole 

Let us investigate the Poynting flux emanating from a rotating quadrupolar magnetic field. The close analytical expression is 
given in terms of spherical Hankel functions such that 


r 4 tt 2 r,2 ■ 2 

Lquad = 5 — cB R sin xi 
3 fM) 


(49a) 


cos 2 \2 


sm X 2 


R 


\h^\k R)\ 2 
1 


H-2“ 

r? 


R 2 


64 


3 

+ T ' 


35 \d r {r 5 \d r (r hfi(k r))\ 


+ 


H— t t: 


\h ( 2 ) { 2 kR )\ 2 r L V 7 \dr(rh { fi( 2 kr))\ 2 R 


In order to get more tractable expressions, we compute the point quadrupole losses given by 


t point _ 128 71 

quad ~ 135 fio c 5 


fi 6 B 2 R 8 sin 2 xi (cos 2 X 2 + 10 sin 2 X 2 ) 


(50) 


It depends on sin 2 xi, which is reminiscent of the dipole field, but now it also depends on the second angle X 2 through a more 
complex trigonometric variation. The difference in spin-down luminosity between the configuration X 2 = 0° and X 2 = 90° 
reaches a factor ten. The exact explicit expression for the true quadrupolar losses being to lengthy, we give the second order 
approximation in a which should be enough in any cases for pulsars. Compared to the point quadrupole, it is given by 

Lquad _ 1 , a 2 (301 - 339 cos 2 x 2 ) 


j point 
quad 


1 + 


+ o(a ) 


(51) 


24 (9 cos 2 x 2 - 11) 

However, the exact analytical variations of L qua d are depicted in Fig.[l] We separated the contributions from the mode m = 1 
(X 2 = 0°) and the mode m = 2 (x 2 = 90°) and assume xi = 90°. Both contributions decrease monotonically with a. Note 
that each luminosity is normalized to its value for a = 0. The largest fluctuations in spin-down are less than 20%. 

To the same order of accuracy in a, the braking index is 


a 2 19 cos 2 X 2 ^ 320 sin 2 X 2 , 2 s 
5 + ) 


_ ( 52 ) 

COS 2 X 2 + 10 sin 2 X 2 “ v “ ' 

It does not depend on xi as it did not depend on x f° r the dipole. It is always less than five for the mode m = 2, decreasing 
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monotonically as shown in Fig. [2] On the contrary, the braking index for the mode m — 1 increases slightly up to n = 5.25 
for a ~ 0.5 before decreasing sharply down ton = 3.17 for a = 1. 


4.3 Hexapole 


The Poynting flux emanating from a rotating hexapolar magnetic field is given by a close analytical expression in terms 
of spherical Hankel functions, see equation m- We do not write down its explicit expression because it is too lengthy. 
Nevertheless, the point hexapole losses can be cast in short form as 


r point _ ® 

hexa ~ 4725 no c 7 


/~v8 rj2 10 • 2 

it B R sin xi 


(53a) 


(29 cos 2 X 2 + sin 2 %2 (1664 cos 2 X 3 + 15309 sin 2 x 3 )) 

It depends on sin 2 xi, but also on the second and third angles (X2,X3) through a more complex trigonometric variation. 
The difference in spin-down luminosity between the configuration ( X2,X3 ) = {(0°, 0°), (90°, 0°), (90°, 90°)} is in the ratio 
29 : 1664 : 15309. These configurations are respectively associated to the mode m = 1:2:3. 

The true hexapolar spin-down is approximately given to second order in a by 

Thexa _ , 7a 2 ( 2 cos 2 X 2 + sin 2 X 2 (60073 — 58025 cos ( 2 x 3 ))) 


j point 
■^hexa 


l - 


+ o(a ) 


(54) 


15 (58 cos 2 X 2 + sin 2 x '2 (16973 — 13645 cos ( 2 x 3 ))) 

The exact variations are shown in Fig. [I] For the mode m = 1, the spin-down slowly increases with a up to 1.04 whereas the 
modes m = 2 and m = 3 decrease with a down to respectively 0.26 and 0.13. 

The corresponding braking index is 

14 a 2 cos 2 \2 + sin 2 X 2 (1024 cos 2 X 3 + 59049 sin 2 X 3 ) 


7- 


+ o(o 2 ) 


(55a) 


15 29 cos 2 X 2 + sin 2 X 2 (1664 cos 2 X3 + 15309 sin 2 xa) 

As shown in Fig. [2] its is less than seven for the modes m = 2 and m = 3, monotonically decreasing to n = 3.77 and n = 3.78 
respectively. For the mode m — 1, the braking index is almost constant, n ss 7, showing a slight increase up to n = 7.13. 


4.4 Octopole 


Eventually, the Poynting flux from a rotating octopolar magnetic field is studied. The point octopole losses are 

16 ,,o 


t point _ 

octo “ 297675 no c 9 


ST 


I T~) 2 r-)12 ■ 2 

B R sin xi 


(56a) 

(56b) 


(23 cos X 2 + sin x '2 (5632 cos X3 + sin X3 (133407 cos x '4 + 11179648 sin X 4 ))) 

The dependence on sin 2 xi is also present, and we retrieve also the second, third and fourth angles (x'2 5 X3 j X&) describing the oc¬ 
topole configuration. The difference in spin-down luminosity between the configuration (x 2 , X3i X 4 ) = {(0°, 0°, 0°), (90°, 0°, 0°), (90°, 90°, 0°), (9 
is in the ratio 23 : 5632 : 133407 : 11179648 (which are for the mode m — 1 : 2 : 3 : 4 respectively). The mode m = 4 has the 
maximal spin-down rate, several orders of magnitude larger than the mode m = 1. Indeed, in all multipole fields, the largest 
spin-down luminosity is attained for l = m. 

The true octopolar spin-down losses to second order in a is given by 

Locto _ , a 2 (68 cos 2 X '2 + 2 sin 2 X 2 (38912 cos 2 X3 + 27sin 2 X3 (1097419 — 999733 cos ( 2 x 4 )))) 


1 - 


£ 0 °“ ~ " 21 (46cos 2 X2 +sin 2 X 2 (11264 cos 2 x 3 + 9 sin 2 X3 (145895 - 116249 cos (2x 4 )))) 


+ o(a ) 


(57) 


The exact variations are shown in Fig. |T] For the mode m = 1, the spin-down significantly increases with a up to 1.24 for 
a Rs 0.7 and then decreases down to 0.95. All other modes m = 2, m = 3 and m = 4 decrease with a down to respectively 
0.59, 0.33 and 0.07. This is a general trend. 

The corresponding braking index is 


:9- 


4a 2 17 cos 2 X 2 + sin 2 X 2 (19456 cos 2 X 3 + sin 2 x '3 (1318761 cos 2 X4 + 28311552 sin 2 x' 4 )) 


+ o(a 2 ) 


(58a 


21 23 cos 2 X 2 + sin 2 X 2 (5632 cos 2 xs + sin 2 X 3 (133407 cos 2 X4 + 11179648 sin 2 X 4)) 

Its exact variations are shown in Fig. [5] The braking index always decreases, whatever the mode m. The final value for a = 1 
can be significantly less than nine. Indeed, for m = 4 the lowest value is n = 4.41 and for m = 3 it is n = 4.34. 

In this paragraph, we demonstrated that the finite size of the magnetic multipole can drastically alter the point multipole 
picture. This reflects in the Poynting flux and therefore also in the derived braking index showing sometimes strong decrease 
with increasing rotation rate. These effects are rooted in the boundary conditions on the stellar surface. We must conclude 
that the electrodynamics on the neutron star surface affects the long term evolution of its rotation. As a matter of fact, 
on non spherical boundary such as those obtained from the electrosphere, results would deviate from the pure point source 
significantly. Explaining the behaviour of intermittent pulsars certainly requires a better and more quantitative treatment of 
those boundary surfaces which should be included with caution. 
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Table 2. Spin-down luminosity for the point multipole source C(l,m ), normalized to its value for m = 1. 


l/m 

1 

2 

3 

4 
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Table 3 . Magnetic multipole strength in Tesla according to the pulsar spin-down luminosity. The fiducial parameters are P = 1 s and 
P = 10 -15 . 

4.5 Implications for the magnetic field strength estimates 

Taking into account the above spin-down luminosities for the multipole fields, we conclude that the magneto-dipole losses 
often used to estimate the magnetic field strength at the neutron star surface is irrelevant if small scale structures exist in 
the crust like the ones induced by high multipole components. The dipole formula only gives a good estimate of the magnetic 
field strength in the far zone, far outside the light-cylinder, but will not furnish a reliable estimate close to the neutron star 
surface where multipolar components are suspected to be dominant. 

The luminosities of the multipoles separated into their components of order (Z, m) are summarized in Table [2] The 
luminosity is normalized according to the mode m = 1. The dependence on the physical parameters is given in 

the last column. In order to get an idea of the misleading field strengths obtained by including only the dipole, we estimate 
the magnetic field on the surface by equating the multipole spin down of order (Z, m) to the rotational luminosity given by 
E = I where I is the moment of inertia of the star. A typical value used here is I = 10 38 kg m 2 . The fiducial pulsar 
parameters used are P = Is and P = 1CP 15 . Results are summarized in Table [3] The dipole is the most efficient radiator in 
the sense that it needs the lowest magnetic field to achieve a certain amount of spin-down luminosity. The next most efficient 
radiator is the quadrupole with l = m = 2. For a given multipole of order Z, the sectorial mode l = m is the most powerful 
radiator. While a dipole needs a magnetic field of the order B « 10 8 T to account for a given spin-down, a quadrupole needs 
B « 10 11 T, a hexapole B ~ 10 14 T and a octopole B « 10 18 T. High multipoles can have magnetic fields above the critical 
field of B q — 4.4 x 10 9 T without contradicting the requirement of dominant dipole losses. Even if multipole spin-down can 
be negligible compared to the dipole, its multipolar magnetic field can be much larger than the dipole field. The upper limits 
for multipole fields Bi tTn seems not to be very restrictive for normal pulsars. 

For millisecond pulsars the situation is different. Having typical period and period derivative parameters of P = 10~ 3 s 
and P = 10 -18 respectively, the constraints on multipoles are tighter. Results are summarized in Table [4] The magnetic field 
is always less than the critical field B q , at least for Z ^ 4. 
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Table 4. Magnetic multipole strength in Tesla according to the pulsar spin-down luminosity. The fiducial parameters are those for 
millisecond pulsars with P = 10 —3 s and P = 10 -18 . 
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log(x) 

Figure 3. Ratio of quadrupole to dipole Poynting flux (blue line), intensity of dipolar field (red line), intensity of quadrupolar field 
(green line) with respect to x for pulsars with a period of one second. Solid lines correspond to X = 1 (m=l) and dashed lines to X = 10 
(m=2). Note the log scales on both axis and that magnetic field strengths are given in Tesla. The black dotted vertical line depicts the 
equality in field intensities Rdip = B quax j. The black dotted horizontal line depicts the equality in Poynting fluxes L<jip = T qU ad- 


The magnetodipole losses formula should be avoided to guess the magnetic field strength at the surface of the neutron 
star. The magnetic field coming from the multipole components are completely missed, although they can be as large or even 
larger than the dipole one, nevertheless without contributing much to the total spin-down power. Observations of the pulsar 
spin-down alone is not reliable to estimate the magnetic field strength at its surface. 

We go on with a discussion about the geometrical effects of multipole fields. Because it becomes very cumbersome to 
investigate in detail the influence of each multipole, having an increasing number of free parameters with increasing l, we 
restrict ourself to the dipole plus quadrupole moments. 

An illustrative example is given by a comparison between the magneto-dipole and magneto-quadrupole losses. We compare 
specifically both m = 1 modes or the m = 1 dipole with the m = 2 quadrupole knowing from the above discussion that 

8tt Bl ip R 6 
3 


Bdip — —— ■ 


Mo c° 


(59a) 


Bquad — 


128 7t x B 2 uad n 6 R 8 


(59b) 


135 mo c ° 

where A' = 1 for m = 1 and X = 10 for m = 2, see table[2] Introducing the ratio between the magnetic quadrupole field B qua d 
and the magnetic dipole field B d i p by 

B qua d 


Bdip 

the associated Poynting flux is a quadratic function of x such that 

Bquad _ 16 X 


R 


(60) 


(61) 


Bdip 45 r£ 

Assuming that the relevant parameters (B, P) are known for each pulsar, we can deduce the dipolar part by fixing x and get 


Bdip — 


3 mo c 3 E 




(62) 


where the spindown is E = 47t 2 I P P~ 3 and / the neutron star moment of inertia taken to be I = 10 38 kgm 2 . We show the 
variations of Bdi p , B qua d and B qU ad/Bdi p for typical pulsar parameters corresponding to a one second period in Fig. [3]and to 
a millisecond period in Fig. [4] Inspecting Fig. [3] for normal pulsars, we recognize that a dominant dipolar Poynting flux does 
not imply a dominant dipolar magnetic field. Indeed for x = 1 the quadrupolar flux is still six to seven orders of magnitude 
less than the dipole flux. They become comparable only for x = 10 3 . If we assume that the Poynting flux remains mainly 
dipolar, we are unable to put severe constraints on the magnitude of multipolar fields such as the quadrupolar components 
(m=l,2). The magneto-dipole losses are useless to get drastic upper limits for multipole components. This conclusion is even 
more severe for magnetars possessing periods longer than one second. 

The situation changes for millisecond pulsars, Fig. [4] For x = 1, dipolar and quadrupolar Poynting fluxes are comparable 
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log(x) 

Figure 4. Ratio of quadrupolc to dipole Poynting flux (blue line), intensity of dipolar field (red line), intensity of quadrupolar field 
(green line) with respect to x for pulsars for millisecond pulsars. Solid lines correspond to X = 1 (m=l) and dashed lines to X = 10 
(m=2). Note the log scales on both axis and that magnetic field strengths are given in Tesla. The black dotted vertical line depicts the 
equality in field intensities f?dip = Rjuad- The black dotted horizontal line depicts the equality in Poynting fluxes L<jip = T qua( j. 


the latter becoming dominant above x = 1. Thus if we assume that the dipole spindown remains dominant, we get stringent 
upper limits for at least the quadrupolar moment, being B qua a < -Bdip- 


5 POLAR CAP GEOMETRY 

It is believed that the geometry of the polar caps determines the radio pulse profiles as the coherent emission is assumed 
to be produced close to the neutron star surface, at heights comparable to the neutron star radius. This region is therefore 
privileged to look for possible significant multipolar components in the magnetic field because a multipole of order l behaves 
like r~ (I+2 i near the origin. 

In order to characterize the quadrupole with respect to the dipole, we introduce two weights wa and w q such that the 
total magnetic field is equal to a weighted linear combination 

B — Wd Bdip "F 1V q Rquad- (63) 

We use the retarded point multipole expressions given in the appendix m Nevertheless, we assume that the neutron star 
radius is equal to R/r l = 0.1 to compute the magnetic field lines. 

To get an idea of the field line topology of the quadrupole compared to the dipole, we show some magnetic field lines 
contained in the equatorial plane for the orthogonal point dipole (x = 90°), fig. [5j and the orthogonal point quadrupole 
((XiiXa) = (90°, 90°)), fig. [U The wave like structure is easily recognized from respectively the double and the quadruple 
spiral pattern. Note however that the magnetic field lines approach only asymptotically these spiral patterns that are not 
shown in the figures for clarity. 

In order to quantify how the quadrupolar field affects the number and geometry of the caps in the vicinity of the surface, 
we show several maps locating these polar caps. The physically relevant geometry depends on the relative weights between 
these dipolar and quadrupolar fields. 

To start with, we recall the polar caps for the dipole field as shown in Fig. [3 for four obliquities: x = {0°, 30°, 60°, 90°}. 
The weights are Wd = l,u> q = 0. The special aligned case with x — 0° represents an azimuthally symmetric situation. Thus 
the polar cap location does not vary with the phase of the pulsar. In the misaligned case, two polar caps are always clearly 
seen. They are directly connected to both magnetic poles. 

Next we consider a quadrupole alone. The weights are Wd = 0, w q = 1. Because this field possesses four poles, we would 
expect at most four polar caps depending on the respective weight of each mode m = {0,1, 2}. A typical sample of polar cap 
maps is shown in Fig. [8] The case xt — X 2 = 90° is the most easiest to interpret. It is a mode m = 2 with magnetic poles 
locate in the equatorial plane. Thus two couples of two poles with opposite polarity are present. These are easily identified in 
Fig. [8] Note however that one cap is hardly seen because of our phase origin which is exactly aligned with one pole. Thus one 
is centred at phase ip = 0. If xt is decreased for instance to xi = 30°, two of these poles move towards the rotation axis while 
the two other stay at the equator, see panel (^i = 30 °,^2 = 90°). Next the case ('xi = 90 °,\2 = 0°) represents the mode 
m = 1 alone. The four caps are still present, but they distribute symmetrically with respect to the equatorial plane. If xi is 
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x/r L 


Figure 5. Geometry of the magnetic field lines in the equatorial plane of an orthogonal point dipole. The dotted circle corresponds to 
the light-cylinder. The field lines approach only asymptotically the double Archimedes spirals (not shown). 



x/tl 


Figure 6. Geometry of the magnetic field lines in the equatorial plane of an orthogonal point quadrupole. The dotted circle corresponds 
to the light-cylinder. The field lines approach only asymptotically the quadruple Archimedes spirals (not shown). 


again decreasing, two caps move to the rotation axis whereas the two other move to the equatorial plane until they merge to 
form only one cap seen at any phase <p, see panel (xi = 30°, \2 = 0°). 

Finally, we assume a linear combination of dipole and quadrupole fields as shown in Fig. [9] For concreteness, we choose 
weights such that Wd = l,w q = 1. Any other weights would also be acceptable but it is impossible to scan all relevant 
parameters and to show the results in a paper of a decent size. Certainly an atlas containing hundreds of maps could 
be computed but we do not show it here. We summarize our results by showing three different classes of behaviour. The 
first class starts with a quadrupole mode m = 2 thus with (xi = 90°, X 2 = 90°), see Fig. [8] We add to it a dipole with 
inclination going from an aligned rotator to a perpendicular rotator. We always find four polar caps located anti-symmetrically 
with respect to the equatorial plane. The configuration becomes exactly symmetric only for a perpendicular dipole, see 
(X = 90°, Xi = 90°, X 2 = 90°). In the other extreme case, see (x = 0°,Xi = 90°, X 2 = 90°), the polar caps are at the opposite 
of what would be expected for an aligned dipole alone. Indeed, their location contradicts the usual fact that an almost aligned 
rotator can only emit one pulse per period. With a quadrupole component, the prediction changes drastically. Two pulses per 
period, separated by half a period, would be detected and most importantly interpreted as a perpendicular (dipole) rotator. 
This is a far reaching important result of this paragraph. If multipole fields exist in neutron star, and we could not understand 
why they should not, all the geometrical fitting parameters obtained by interpreting the radio pulses would be inaccurate or 
worst, wrong. 

The second case assumes a quadrupolar symmetric mode m = 0 to which we add a perpendicular dipole as seen in the 
panel (x = 90°, xi = 0°,X2 = 0°). This configuration depicts the opposite example of the previous case, in which a clear 
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phase (ip) 


phase ( ip ) 


Figure 7. Geometry of the polar cap for the retarded point dipole field solely with R/r l = 0.1. The four panels correspond to different 
inclination angles \ = (0°, 30°, 60°, 90°} as indicated in the plots. 


misinterpretation of perpendicular rotator being an aligned rotator is performed. We always only see one pulse per period 
with a duty cycle between 20% and 50%. This behaviour is usually explained by an almost aligned rotator with \ close to 
zero. To the contrary, here, the dipole is perpendicular but the quadrupole is aligned. 

The third and last case corresponds to the quadrupole mode m = 1. Thus we add a m = 1 dipole to a m = 1 
quadrupole. We always observe four polar caps. For too small obliquity x 45°, on polar cap disappears, only three remain, 
see panel (x = 0°,Xi = 90°, = 0°). For sufficiently large obliquity \ > 45°, polar caps located in the same hemisphere 

become identical. 

Before dealing with the high-energy emission properties, we stress the main important result of this paragraph. If pulsar 
magnetospheres are only described with a misaligned dipolar magnetic field, interpretation of the radio pulse profile becomes 
extremely misleading for some special cases. Taking into account quadrupole or higher order multipole fields would destroy 
the fitting obtained by a dipole solely and the value of the parameters would be questionable. 


6 HIGH-ENERGY EMISSION PHASE DIAGRAM 


The multipolar magnetic field impacts the emission in the whole electromagnetic spectrum, not only radio as shown in the 
previous section but also the high energy counterpart. This is especially true for millisecond pulsars, those having a large R/r l 
ratio. In order to show the complexity of light curves emanating from such magnetic configurations, we plot several phase 
diagrams including dipolar and quadrupolar fields posses sing s i milar relative strengths (wa = w q = 1). In order to compute 
the light-curves we used a slot gap model as described in IPetri (2015). 

The phase diagram for the dipole is shown in Fig. [10] The caustic effect is seen in red-blue-black colour, corresponding 
to the highest intensity. The non emitting phases remain in white. They are larger than the related polar caps because the 
magnetic field lines are not directed radi ally but their opening angle is larger than the cone subtended by the polar cap by 
approximately 50% (Ganeadhara fe Gupta 200 1!'). 

For the quadrupole alone, we recognize an m = 1 structure for (xi = 90°, = 0°) with two main peaks whereas for a 

m = 2 structure with (xi = 90°, = 90°) we observe four peaks, see Fig. [TT] Two other intermediate cases are also shown 

corresponding to the same configuration as in Fig. [8] Note that in general the north-south symmetry is broken by the presence 
of a quadrupole. The case (xi = 30°, X 2 = 0°) possesses a significant bridge emission between both pulses around (/ = 90°. 
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phase (tp) 


phase (p) 


Figure 8. Geometry of the polar cap for the quadrupole field solely with R/rt, = 0.1. The four panels correspond to different configu¬ 
rations depicted by the couple of angles (xi>X2) as indicated in the plots. 


The most interesting case uses a linear combination of dipole and quadrupole as done in equation (16311 . The variety of light- 
curves is shown in Fig. [12] where the caustics are identified in red-blue-black colour. All the possible light-curves are present, 
four pulses, three pulse, two pulses, one pulse, and sometimes with bridge emission. Another sample is shown in Fig. 1131 A 
degeneracy between the angles (XiXIjXz) shows that different combinations leads to the same phase diagram, modulo some 
translation in phase or reflection with respect to the equatorial plane. This is the case for (x = 90°, xi = 180°, X 2 = 0°) 
which is very similar to (x = 90°, xi = 0°,X2 = 0°) apart from a translation of half a period. The same applies for 
(X = 30°, xi = 180°, X 2 = 0°) and (x = 30°, xi = 0°,X2 = 0°) which are reflection symmetric to each other with respect to 
£ = 90° and translated by half a period. This expectation is immediately seen from the chosen expressions for the constants 
<li,m possessing themselves some (anti-)symmetry properties. 

Symmetry considerations can significantly decrease the number of relevant plots to show in a an atlas. 


7 CONCLUSION 

It is hard to believe that a neutron star only possess a dipolar magnetic field. Although the rotating dipole is the dominant 
component in the far wave zone, higher order multipoles must be present close to its surface and in its magnetosphere. We 
solve analytically the time harmonic Maxwell equations in vacuum and the corresponding expression for the Poynting flux 
taking into account the finite size of the star. We derived useful analytical expressions for the Poynting flux, the braking 
index and the retarded point multipole fields. We briefly discussed the important implications for magnetic field estimates 
at the neutron star surface, the radio pulse and high energy emission properties. We showed that already in the vacuum 
case, expectations and predictions relative to a pure dipole are drastically diverging. For instance, a double peaked light 
curve separated by approximately half a period can be explained by a geometric configuration deviating significantly from an 
orthogonal rotator if we accept the presence of a quadrupolar field of the same intensity as the dipole. 

We plan to look more deeply into the phase resolved radio polarization characteristics to constrain the field line geometry 
close to the star. Including our picture in a force-free pulsar magnetosphere would also alleviate the conclusions drawn so far 
in the literature about pulsar electrodynamics in a cold massless plasma. 
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Figure 9. Geometry of the polar cap for a linear combination of the dipole and quadrupole fields with = 0.1 and weights 

w d = w q — 1- Each panel corresponds to a different configuration depicted by the angles (x- Xl ■ X 2 )- Their precise values are indicated 
in the plots. 
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Figure 10. Phase diagram for the retarded point dipole rotator with R/rj^ = 0.1. The four panels correspond to different inclination 
angles x = {0°, 30°, 60°, 90°} as indicated in the plots. 
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Figure 11. Phase diagram for the retarded point quadrupole rotator with R/r l = 0.1 and inclination angles (xi 5 X 2 )- 
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Figure 12. Phase diagram for the retarded point dipole and quadrupole and the inclination angles (x, Xi ■ X 2 ) for weights = w q = 1. 
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APPENDIX A: COMPARISON WITH PREVIOUS WORK 


As a check of our formalism, we compare our approach to the more cumbersome technique used bv lBonazzola et al.l (120151 1. 
For the axisymmetric multipole magnetic field, if we substitute their definition of the constant Bj 3 by 


r i+2 b o = _ a L 
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1(1 + 1 ) (21 + 1 ) 
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we retrieve the magnetic field given in their equation (16). For the electric part we get 
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This is the opposite of the electric field given in equation (21) of Bonazzola et al. ( 2015 1 (apart from the monopolar term 
that we did not include in our work). Indeed their E-s component is not consistent with their boundary condition (19) and 
recurrence formula (20). It should be the opposite sign as given above. 

For the asymmetric case, to retrieve the expansion into transverse electric (TE) a nd transverse magnetic (TM) modes, 
we should make the following identification between our notations and the one used by Bonazzola et al. ( 20151 1 


Al,m — —km Cl+m +1(1+1) 


E0 Ai m — —km a+m +1(1+1) 


21 + 1 (1 — m)\ 
4n (1 + m)\ 

21 + 1 (1 — m)\ 


(A4a) 

(A4b) 


4n (1 + m)\ 

Therefore the TE mode is associated in our picture to the /, s m part of expansion in equation (HO)-® and the TM mode is 
associated to the f im part of the same expansion. It is straightforward t o che ck that we indeed get the magnetic and electric 
field components as reported in equations(26)-(29) of Bonazzola et al. ( 2015ll . The last step concerns the imposition of the 
boundary conditions on the neutron star surface. Adopting their definition Bf. lm of the internal field we have 


p< (R \- /,(, | n / 2 1 + 1 (l~m)\ fi B m(R) 

B r;lm(R) ~ " V* Q + V ~ 


Using our definition of /; S m (R) from eq. (1101) and the identification of the TE modes in ea. (IA4all we get 


m- 


km P 


(km R) 

which is nothing but equation (30) of Bonazzola et al. ( 2015l l. For the TM modes, using the identification eq. (IA4bl) and our 
boundary conditions eq. m we get 
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These are identical respectively to equation (32) and equation (33) of Bonazzola et al. ( 2015l l. Consequently, we demonstrated 
that our formalism employing vector spherical harmonics gives exactly the same results as a decomposition into TE and TM 
modes. 


APPENDIX B: RETARDED POINT MULTIPOLES 
B1 Retarded point dipole 

Boundary conditions on the neutron star reflect in the vacuum solution, even very far away in the wave zone, through additional 
1 = 2 terms in the electromagnetic field. Their effect is dominant for extremely rapidly rotating stars, a = R/r l < 1. As 
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neutron stars are observed to rotate at most with a period of 1.6 ms, the largest ratio of R/rr, corresponds roughly to 0.1. 
The Deutsch solution can then be expanded in a rapidly converging series of R/r l. For practical purposes, it is a very good 
approximation, even to the lowest order. To get this approximation, we compute the solution in the limit of vanishing stellar 
radius by setting R —>• 0 in the Deutsch solution given for instance in full length in IPetril (I2012h . Nevertheless we assume a 
constant magnetic moment p = B R 3 in this limit. We obtain the following expressions for the component of the magnetic 
field in spherical coordinates 


B r = 


B-a — 


2 B R 


[cos x cos i? + sin x sin d (cos ip + kr sin ip)] 


BR: 


B v = — ■ 


[cos x sin d + sin x cos d {( k 2 r 2 — 1) cos ip — kr sin ip)] 
sin x [k r cos ip + ( k 2 r 2 — 1) sin ip] , 


(Bla) 

(Bib) 

(Blc) 


with the instantaneous phase given by -i p = kr — Qt + ip. We draw the attention of the reader not to confuse the above retarded 
point dipole solution (ED with the exact Deutsch field given in equation m- In the point dipole solution, the quadrupole 
terms induced by the corotating quadrupole electric field inside the star are missing because the boundary conditions are 
irrelevant for a point dipole (so no cos 2$ term for instance). In general the retarded point multipole is not exactly the same as 
the finte size multipole field, the latter taking into account boundary conditions on a sphere of finite size. In a point multipole 
description, there is no need to impose boundary conditions, therefore they are no higher multipolar electric field generated 
by the rotating magnetic field because no charge redistribution in a perfectly conducting sphere. Nevertheless the point dipole 
approximation is often used to model pulsed emission from neutron star magnetospheres. We found it useful to give simple 
analytical expressions for any retarded point multipole in the limit of vanishing r adius of the spher e. It is computationally 
more efficient without losing accura cy (to second order in R/rp). See for instance iDvks et al.l (120041 1 for the retarded point 
dipole in cartesian coordinates and iMelrose fe Yuen ( ,20141 1 on how to derive this simple point dipole solution with vector 
algebra. These components are used to compute magnetic field lines in pulsar dipolar magnetospheres to get the shape of 
the polar caps, see [JS]and from which pulsed emission is supposed to come from, see m The same procedure as before can 
be reproduced for higher order multipoles. We show three more examples concerning the quadrupole, the hexapole and the 
octopole in the next paragraphs. 


B2 Retarded point quadrupole 


Expanding the exact vacuum solution to lowest order while keeping the magnetic quadrupole moment constant, we can derive 
the components of the magnetic field. In this case, we show separately the contribution of each mode m. Note here again that 
the expressions given below apply only to a point quadrupole, without specifying boundary conditions on the neutron star 
because its size is zero in this limit. The axisymmetric mode m = 0 has components with weight cos;\;i such that 
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The first asymmetric mode, m = 1, has components with weight sinxi COSX2 such that 
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— \ — cos d(3(5 — 4 k 2 r) sin ip + fcr(7fc r 2 — 15) cos ip) 

Q W 25 


The last mode m = 2 has components with weight sin xi sin \2 and are given by 

- sin 2 d((4fc 2 r 2 — 3) cos 2ip — 6 kr sin 2 ip) 


B r = 

BR 4 


r 4 


B$ = 

BR 4 

1 

r 4 

3 

Bcp = 

BR 4 

1 

r 4 

3 


— \j — sin 2d(2fcr(3 — 2k r 2 ) sin 2ip + (3 — 6 k 2 r) cos 2ip) 
^\/ sind(3(2fc 2 r 2 — 1) sin 2 ip + 2fcr()3 — 2k 2 r 2 ) cos 2ip) 


(B2a) 

(B2b) 

(B2c) 

(B3a) 

(B3b) 

(B3c) 

(B4a) 

(B4b) 

(B4c) 
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Multipolar electromagnetic fields 27 

These expressions together with the retarded point dipole are used to compute the polar cap shape as well as the high-energy 
pulsed emission in the multipole-caustic model described in ||5]and |J 6 ] 


B3 Retarded point hexapole 


By expansion of the exact vacuum solution to lowest order, keeping the magnetic hexapole moment constant, we show 
separately the component of each mode m. Note here again that the expressions given below apply only to a point hexapole, 
without specifying boundary conditions on the neutron star because its size is zero in this limit. The axisymmetric mode 
m = 0 has components with weight cos Xi such that 


B r = ——t— 771 / 77(3 cos i? + 5 cos 3fl) 
r a 2 \ 2 

B$ = — —pf ^ y^"(sin i? + 5 sin 3i?) 

B v = 0 . 


(B5a) 

(B5b) 

(B5c) 


The first asymmetric mode m = 1 has components with weight sin xi cos X 2 cos ip given by 


B r = — 
B$ = — 


B R b 1_ 
r 5 10 

B R 5 cos i?(37fc 4 r 4 + 183fcV + 35 cos 2i?(fc 4 r 4 - 21fcV + 45) - 735) 


%/21sini?(5cos2i? + 3) (2 k 2 r 2 — 5) 


40^21 


Btp = 


B R 5 kr (5 cos 2i?(46fc 2 r 2 - 105) + 42fc 2 r 2 - 315) 


40^21 


and components with weight sincos X 2 sin ip given by 


B r = —^-\/-^fcrsini?(5cos2i? + 3 )(k 2 r 2 — 15) 

r° 10 V 3 


B$ = 


Bco = 


B R b kr cos i?(105 cos 2 i ?(2 fcV - 15) + 62fc 2 r 2 + 735) 
r 5 40\/2T 

BR 5 5cos2i?(13fc 4 r 4 - 81ft V + 105) + 7(fc 4 r 4 - 21 ft V + 45) 


40\/2l 

The mode m = 2 has components with weight sin xi sin X 2 cos X '3 cos 2ip 

BR 5 [42 . 


B r = —4— \l -T- sin 2 i9 cos i?( 8 k 2 r — 5) 

B R 5 sin i?(7 cos 2i?(16fc 4 r 4 - 84ftV + 45) + 3(48fc 4 r 4 - 92ft 2 r 2 + 35)) 


Btf = 


B v = — 


4\/2l0 


B R 5 fcr sin 2i?(76fc 2 r 2 — 105) 
r 5 fi210 


and components with weight sin \i sin X 2 cos X 3 sin 2 ip 


B r = 


B R b /14 


Btf = — 


kr sin i? sin 2i?(4fcr 2 — 15) 

BR 5 kr sin i?(21 cos 2i?(8fc 2 r 2 - 15) + 136 k 2 r 2 - 105) 


2V2W 


B R 


18 \l -pp:k 2 r 2 sin2 1 ?. 


Finally, the last mode m = 3 has components with weight sin^i sinx 2 sinx 3 cos3i p is 

B R 5 f7 . 


Br =- 


B$ = — 


j sin 3 i?(18fe 2 r 2 — 5) 
o 


r 5 4 V 5 

_ BR 5 9 fir 

”'*= — r 5 4 V 5 


— sin 2 i?cos i?(9fc 4 r 4 — 21fcr 2 + 5) 


— \ —kr sin 2 i?( 6 k 2 r — 5) 


(B 6 a) 

(B 6 b) 

(B 6 c) 

(B7a) 

(B7b) 

(B7c) 

(B 8 a) 

(B 8 b) 

(B 8 c) 

(B9a) 

(B9b) 

(B9c) 

(BlOa) 

(BlOb) 

(BlOc) 
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28 J. Petri 


and components with weight sinxi sinx 2 sinx 3 sin 3*0 

BR 5 l~7 . 3 q / q ,2 2 r\ 

B r = -=— o\l — hr sin r — 5) 

r 5 V 5 

B R 5 9 [ 7 1 . 2 Q q ^,2 2 

= —r- \ — At sin vcosmbk r — 5) 

v b 4 V 5 

B R 5 3 /~7 . 2 Q/ n7 4 4 o-, 7 2 2 , 

Ben = —--- \ — sin $(9« r — 21/c r + 5) 

f 3 4 V o 


(Blla) 

(Bllb) 

(Bile) 


B4 The retarded point octopole 


Performing the lowest order expansion of this octopole, we get the retarded point octopole which is splitting into its various 
azimuthal modes m. Note here again that the expressions given below apply only to a point octopole, without specifying 
boundary conditions on the neutron star because its size is zero in this limit. Explicitly, the mode m = 0 has components 
with weight cos yi and 


B R ° 1 / I e: 

B r = -— — \ —(20 cos 2$ + 35 cos 40 + 9) 

r e 32 V 2 v ’ 

B R 6 1 /15, . _ _ . 

B$ = -— — W — (2 sm 20 + 7 sm 40) 

= 0 . 


The mode m = 1 has components with weight sin \i cos \2 cos ip 
D _ B7 ? 6 (2 sin 20 + 7 sin 40)(fc 4 r 4 — 45k 2 r 2 + 105) 

U271 

BR 6 4cos20(19fc 4 r 4 — 73k 2 r 2 + 63) + 21(fc 2 r 2 (2fc 2 r 2 — 5) + cos40(2fc 4 r 4 — 39k 2 r 2 + 84)) 
r 6 336^3 

B R 6 kr cos 0(—59fc 4 r 4 + 795fc 2 r 2 + 21 cos 20(7fc 4 r 4 — 135fc 2 r 2 + 315) + 945) 

V _ r 6 1260%/3 


and components with weight sincos \2 sin ip 


D BR 6 5kr (2sin20 + 7sin40) {2k 2 r 2 - 21) 

^ lWf 

B R 6 kr 

BrQ = ---— (4cos2 , $(46/c 4 r 4 — 780k 2 r 2 + 945)+ 

r 6 5040^3 


21(5fc 2 r 2 (fc 2 r 2 — 15) + 3cos40(fc 4 r 4 — 55k 2 r 2 + 420))) 

_ BR 6 cos0(—23fc 4 r 4 + 32k 2 r 2 + 21 cos20(3 k 4 r 4 - 16fc 2 r 2 + 21) + 63) 
~~ r 6 84^3 


The mode m — 2 has components with weight sin xi sin \2 cos X 3 cos 2 ip 


B r = 

Btf = 

B v = 


BR 6 sin 2 0(7cos20 + 5)(16fc 4 r 4 - 180fc 2 r 2 + 105) 
r 6 28p6 

B R 6 sin 20(88 k 4 r 4 + 47fc 2 r 2 + 21 cos 20(8fc 4 r 4 - 39 k 2 r 2 + 21) - 63) 
r 6 42^6 

B R 6 kr sin 0(21 cos 20(32fc 4 r 4 - 240fc 2 r 2 + 315) + 5(64fc 4 r 4 - 600fc 2 r 2 + 945)) 
r 6 315^6 


and components with weight sin yi sin \2 cos \3 sin 2 ip 


B r 




B v 


B R 6 5 kr sin 2 0(7cos 20 + 5)(8 k 2 r 2 — 21) 

14^6 

B R 6 kr sin 20(21 cos 20(4fc 4 r 4 - 55k 2 r 2 + 105) - 5(37 k 2 r 2 + 63)) 

105%/6 

B R 6 sin 0(176fc 4 r 4 - 620fc 2 r 2 + 21 cos 20(16fc 4 r 4 - 44fc 2 r 2 + 21) + 315) 
r 6 42^6 


(B12a) 

(B12b) 

(B12c) 

(B13a) 

(B13b) 

(B13c) 

(B14a) 

(B14b) 

(B14c) 

(B15a) 

(B15b) 

(B15c) 

(B16a) 

(B16b) 

(B16c) 
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Multipolar electromagnetic fields 29 


The mode m = 3 has components with weight sin xi sin \2 sin X 3 cos X 4 cos 
= 1 \j sin 3 •& cos #(27 k 4 r 4 — 135k 2 r 2 + 35) 


B 

B$ = — 


r ~ 2 V 7 ' 

B R 6 sin 2 #(144fc 4 r 4 - 193fc 2 r 2 + 3 cos 2#(54fcV - 117fc 2 r 2 + 28) + 42) 

4%/2T 


r-6 


B v = — - 


= -sin 2 $ cos #(87fc 4 r 4 — 415fc 2 r 2 + 315) 


r u ll) v i 

and components with weight sin xi sin X 2 sin X 3 cos X 4 sin 3tp 

D E > 6 ic AT 


=- 


B R 6 15 /If, . 3 q ^ (al 2 2 
—---W-At sin VQ,os r u(pk r — 7) 


B# = -U/ffcr sin 2 #(93fcV - 335fcV + cos 2#(81fcV - 495fcV + 420) + 21( 

r b 20 V 7 

BR 6 sin 2 i? cos #(153A; 4 r 4 - 272k 2 r 2 + 63) 

2fi21 


J-f ± L Dili U V^V_»0 V I / Zj | Zj I4j I I KJtJ I 

The mode m = 4 has components with weight sin xi sin X 2 sin X 3 sin X4 


cos Aip 


B r 

B$ 


BR 6 sin 4 #(256fc 4 r 4 - 720fcV + 105) 
r 6 4^42 

BR 6 sin 3 i?cos$(128fc 4 r 4 - 156/cV + 21) 
r 6 y[42 


BRi 


■ — \j-pfikr sin 3 $(64fc 4 r 4 — 220k 2 r 2 + 105) 
o V J- 


and components with weight sinxi sinx 2 sinx3 sinx4 sin 4^ 

BR 6 5£x sin 4 #(32fc 2 r 2 — 21) 


B r = 


V42 


Btp 


—-tt— —\ — kr sin 3 # cos #(64fc 4 r 4 — 220k 2 r 2 + 105) 
r b 5 V 21 

BR 6 sin 3 #(128fc 4 r 4 - 156fcV + 21) 
r 6 -^42 


(B17a) 

(B17b) 

(B17c) 

(B18a) 

(B18b) 

(B18c) 

(B19a) 

(B19b) 

(B19c) 

(B20a) 

(B20b) 

(B20c) 
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